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We apply the SU(2) slave fermion formalism to the Kitaev honeycomb lattice model. We show that both the 
Toric Code phase (the A phase) and the gapless phase of this model (the B phase) can be identified with p-wave 
superconducting phases of the slave fermions, with nodal lines which, respectively, do not or do intersect the 
Fermi surface. The non-Abelian Ising anyon phase is a p + ip superconducting phase which occurs when the B 
phase is subjected to a gap-opening magnetic field. We also discuss the transitions between these phases in this 
language. 



I. INTRODUCTION. 

In Ref. Q] Kitaev introduced the following remarkable 
model of s = 1 /2 spins on a honeycomb lattice 

H = -J x SjS? - J v 2J SjSj — J* SjSj > 

x— links y— links z— links 

(1) 

where the z-links are the vertical links on the honeycomb lat- 
tice, and the x and y links are at angles ±7r/3 from the verti- 
cal. This model is exactly solvable and has a gapped Abelian 
topological phase (the 'A phase') which is equivalent to the 
Toric Code 2 . It also has a gapless phase (the 'B phase') which, 
when subjected to an appropriate time-reversal symmetry- 
breaking perturbation, becomes a gapped non-Abelian topo- 
logical phase supporting Ising anyons. 

This model is one of the rare instances of an exactly solv- 
able model of a quantum magnet which does not order in 
its ground state and, instead, condenses into a topological 
phase. As such, it is a useful testing ground for theoreti- 
cal techniques, such as slave fermion representations, which 
have been applied to approximately solve models of frustrated 
magnets which are not exactly solvable. Applying these tech- 
niques to Eq. [T] can shed light on the physics of this model 
and, conversely, on the applicability of these techniques. 

Kitaev solved the Hamiltonian ([T]l by introducing a 
fermionization of the spins in terms of Majorana fermions. 
By expressing each spin operator as a product of two Ma- 
jorana fermions, the spin model can be described exactly as 
a model of Majorana fermions coupled to a Z2 gauge field. 
In this description the effect of the gauge field is particularly 
transparent: the physical correlators are captured exactly by 
the fermionic band structure, and the gauge field serves only 
to enforce the fact that only gauge-invariant observables (e.g. 
products of spins) are physical. 

In this paper, we apply a different fermionization proce- 
dure, the SU(2) slave fermion formalism. This representation 
requires a different projection to eliminate redundancies in the 
Hilbert space compared to Kitaev's representation in terms of 
Majorana fermions; therefore, it is interesting to see how the 
same low-energy degrees of freedom emerge. In the SU(2) 
slave fermion formalism, the spins are written in terms of stan- 
dard, rather than Majorana, fermionic spinons. The Hamilto- 
nian of Eq. [T]is then expanded about an RVB mean field state. 



We show that this is a stable mean-field theory which captures 
the the physical correlation functions of the exact ground state 
of Eq. [T] We find that the A phase is a p-wave supercon- 
ducting state of the slave fermions. The state is fully gapped 
because the nodes in the order parameter do not intersect the 
Fermi surface. The Majorana fermions of Kitaev's solution 
appear as Bogoliubov-de Gennes quasiparticles of the super- 
conducting state. The B phase is a p-wave superconducting 
state with gapless excitations at the nodal points. These exci- 
tations form a single Dirac fermion. When the order param- 
eter develops an ip component, the Dirac fermion acquires a 
mass, and the system goes into an Ising anyon phase. The 
transition point between the A phase and the gapless B phase 
is an interesting quantum critical point, which we describe in 
terms of superconducting order parameters. 

By studying the theory of fluctuations about the mean-field 
saddle point, we recover the Z 2 gauge field as the unbroken 
gauge symmetry remaining in the superconducting state. This 
situates the ground state of the finely-tuned Hamiltonian ([T]) 
in the broader context of spin liquid 3 10 and superconducting 
phases, and allows us to understand its phase diagram in terms 
of these more familiar phases of matter. 

II. SU(2) SLAVE FERMION FORMULATION 
A. Slave Fermion Mean-Field Hamiltonian 

Our starting point is the representation of the spin operators 
in terms of spinful Dirac fermions, first discussed in Ref. |3] 
We thus write the spin operators Sj, i = x,y, z, as: 

Here, we have introduced the fermion operators /j Q , usually 
called spinons. For two-spin interactions of the form Sf§j, 
one way to treat the resulting Hamiltonian is to use a Hubbard- 
Stratonovich transformation to decouple the 4-fermion inter- 
actions, re-expressing them as interactions between a bosonic 
field $ (which lives on a link in the lattice) and a pair of 
fermion operators on the sites i and j bordering this link. 
One may then study the mean-field solutions which can be ob- 
tained by condensing the bosons. This is often a fruitful way 



2 



to investigate candidate 'spin liquid' ground states, in which 
the spins are strongly correlated but have no spatial order. 

One important caveat in this formulation is that Eq. [2]gives 
a faithful representation of the Hilbert space only in the sub- 
space of fermionic states for which each site is singly occu- 
pied. Thus at each site (i), we must impose the 3 (redundant) 
constraints 



grangian on the .r-links can be written in the form: 



f t f t _ 
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Hit + "-4 = 1 
fitfil = ■ 



(3) 



As explained in Ref. 01 when the Hamiltonian preserves 
SU(2) spin rotation symmetry, the Lagrange multipliers of 
these constraints can be viewed as the temporal component 
of an SU(2) gauge field, leading to a theory of fermions cou- 
pled to a fluctuating gauge field. (The spatial components of 
this gauge field are given by the phases of the fermion kinetic 
terms, which here arise due to condensation of a bosonic field 



- see Appendix B 3 i. Projection would be enforced by inte- 
grating out the gauge fields. In practice, this is typically done 
approximately using perturbation theory in the fermion-gauge 
field coupling 2 ^. 

Thus the decoupling (|2j leads to a description of the spin 
model as a theory of fermions (spinons) coupled to an SU (2) 
gauge field. For the Hamiltonian ([TJ, we will find that the 
spinons are in a superconducting phase, such that this gauge 
symmetry is broken down to Z2, and in particular is fully 
gapped, such that the effect of dynamical gauge-field fluc- 
tuations on the fermion band structure is minimal. We will 
nonetheless find that this gauge theory is a useful tool to un- 
derstand the origin of the various topologically ordered phases 
described in Ref. [TJ 

We begin our analysis with the mean-field description of 
the exact spin-liquid ground state of the Hamiltonian ([TJ. In 
the case of spin-rotationally-invariant Hamiltonians, such as 
the Heisenberg model, the Hamiltonian simplifies consider- 
ably when written in terms of the fermions |2j. In the absence 
of spin-rotational symmetry, as in Eq. [TJ the Hamiltonian is 
more complicated. For instance, the Hamiltonian on x-links 
takes the form. 
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+ *i (firfit + fljji ) + [ frfr - flfii ) + h.c. 

+ e a + fljl) + ;<->, (44 - fyi) + h.c. 

where h.c. is the hermitian conjugate with all spin directions 
reversed. The Lagrangian can be decoupled in a similar man- 
ner on the y- and x-links as well, as detailed in Appendix [B] 

Before proceeding, it will be helpful to pick a unit cell for 
the honeycomb lattice. We will label the two different sites 
with a unit cell by the index i = 1, 2 and different unit cells 

by R = nix + H2(|x + ^y). Then, we denote the fermion 
fields by fnia- Their Fourier transforms are defined by: 
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N 

R 



(5) 



where N is the total number of lattice sites. 

To proceed, we assume that acquire non-zero ex- 

pectation values. We parametrize these expectation values by 
Uj ta , Aij_ a , a =t, as explained in Appendix B2 Unlike 
in the case of Heisenberg interactions, to describe the Kitaev 
model we must condense both hopping and superconducting 
order parameters or else the mean-field equations will not be 
satisfied (except in the special case J x = J y = 0, J z 7^ 0), 
as shown below. (In the Heisenberg case, hopping and d- 
wave superconducting terms can be rotated into each other 
by a gauge transformation. This is not true for the p-wave 
superconducting case considered here.) Because SU(2) spin 
rotation invariance is explicitly broken on each link, the lat- 
ter involve the spin-polarized superconducting terms A-j-, A^. 
Thus, replacing the fields <£>j, 8; by their expectation values, 
we obtain the mean-field Hamiltonian: 



H 



Mq) A ff (q) - 

C(q) -A CT (-q) 

-A*(- q ) -C(-q) 

A*(q) -M-q) 



q,er 

^q = ( /q,l,CT /q,2,<r /-q,l,ff /-q,2 



■0q<T 

(6) 



(4) 



(Here the factor of | in the first line compensates for the fact 
that the expression |6j counts each term in the Hamiltonian 
twice. Alternatively, we could sum over half the Brillouin 
zone.) If we write ip q in components, it has three indices (in 
addition to momentum), ipqi aa , where i = 1, 2 is a sublattice 
index, a =j, i is a spin index, and a = ± is a particle-hole 
index. 

Since we will often be using Pauli matrices to act on these 
1 indices, we will, to avoid confusion, introduce three different 

Si Sj + S { Sj + S i Sj = —yfiafjafjpfip + jfiafiafjpfjP notations for Pauli matrices. We will use a^ ,z for Pauli ma- 
trices acting on spin indices; [i*j V ' z for Pauli matrices acting 
on sublattice indices; and t^ v ' z for Pauli matrices acting on 
particle-hole indices. (Of course, it is precisely the same three 
matrices in all three cases.) 



with similar terms on the y-links, as detailed in Appendix 
[B] (This form is not unique; using the constraints, it can be 
rewritten in different forms which are equivalent in the con- 
straint subspace.) In the Heisenberg model, by contrast, the 
Hamiltonian on each link can be written in the form: 



1 



ft 



1 



■t 



As a result of the more complex form of the Hamiltonian, 
it is necessary to introduce four Hubbard-Stratonovich fields 
to decouple the four-fermi interactions. For example, the La- 
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By requiring self-consistency of the expectation values, we 
can express Uj iCC , A.y. Q in terms of J x ,y,z> as shown in Eq. 
B7 At the saddle point of interest, the relevant parameters 
are: 



*t(9) = 




A t (g) = 




h(q) = 




A+(<z) = 





J, 



J, 



J, 



(7) 



where fi.2 = ± |i are the lattice vectors. 

The band energies and eigenfunctions of Hmf reveal 
the correspondence between this picture and the Majorana 
fermion decoupling of Ref. 1 The mean-field spectrum con- 
sists of 3 flat bands, with energies: 



± J* ( 



± Jy ( 



± 



Jz 



(8) 
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(11) 



where 9„ = 



Jll 



J z , and in all cases 



+ corresponds to the negative-energy solution. The ,- L there- 
fore lie in the 3 fiat bands, and are localized on x,y, and 
z links respectively, and c is the dispersing Majorana mode 
identified by Ref. [JJ 

Hence the saddle point |7]i reproduces exactly the descrip- 
tion of Ref. 1 with the precise mapping between the fermions 
' g)OV s and Kitaev's Majorana fermions given by Eq. 
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only difference is that Ref. 1 does not include the energy of 
the flat bands, so that b x ' y,z enter only in determining the band 
structure of the remaining Majorana mode c. The fermionic 
mean-field energy we obtain per unit cell at half-filling 



1 



(J x + Jy + J Z ) - 



(12) 



and one dispersing band, of energy 



1 



e^q) = ±-\J x e^ h + J y e 1 ^ + J z \ 



(9) 



(Since we have included an explicit factor of 1/2 in our def- 
inition of the spin operators Si, our J x ,y,z are 4 times larger 
than Kitaev's. There is an additional explicit factor of 4 in his 
definition of the spectrum in Eqs. 31 and 32 in Ref. Q] This 
accounts for the factor 16 between our spectra.) The corre- 
sponding eigenvectors are naturally expressed in terms of the 
Majorana fermions 



However, the first term is cancelled by the zero-point energy 
arising from terms of the form 



. , in the Hubbard- 
Stratonovich Hamiltonian, so we are left with precisely the 
same energy as in Kitaev's solution. 

Superficially, we have obtained an 8-band mean-field the- 
ory from a model of spinful fermions on a lattice with a 2-site 
unit cell. Readers might thus justifiably be concerned that we 
have in fact obtained double the degrees of freedom that we 
would have expected. However, we have combined / q i CT and 

f-qicr l nt0 tne same spinor; consequently, we should restrict 
q to half the Brillouin zone to avoid double-counting. 
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(/q4 f- 



-q4 



(10) 



We have used the same labels as Ref. Q]for these operators. 

However, this is not a unique mapping. For instance, we 
could, instead, take c = -(fl + fa), b x = i(fl - /j,), 

b v = /| + fi, b z = — ff). Furthermore, the mean- 
field Hamiltonian has a different expression in terms of these 
operators than in the mean-field theory of Ref. QQ For exam- 
ple, the bilinears b z R1 b z R2 do not commute with the mean-field 
Hamiltonian. The reason for this is that if the spin operators 
are expressed in terms of the /, ps according to Eq. [2] and 
then the /, ps are written in terms of c, b x ,b y ,b z , according 
to Eq. [TO] then we will not obtain the same representation 
as in Ref. OQ Only after the constraints are imposed do the 
operators in Eq. [10] become equivalent to Kitaev's. This is 
explained in more detail in Appendix [A] 

The eigenvectors corresponding to the eigenvalues ([8]) and 
(|9]l are given by: 

a-±(«) = \ {i^Hl x ±bl 2 ) 



B. Slave Fermion Band Structure 

To understand the physics of this model, it is useful to focus 
on the band structure of the down-spin fermions. It suffices to 
consider the case J x = J y = J: 




. . ( q x . V3q v 

+4 cos — sin 

I 2 2 



1/2 



(13) 



This describes a pair of bands which cross at either or 2 
distinct points in the Brillouin zone. Following Ref. UJ we 
will call the former case, which occurs for \J Z \ > 2\J\, the A 
phase. In the A phase, the spectrum is fully gapped. When 
\J Z \ < 2\J\, there are two Majorana cones in the spectrum or, 
equivalently, a single Dirac cone. This is the B phase. Our 
objective here is to understand how this band structure arises 
in the slave fermion superconductor, and use this analogy to 
understand the transitions between these phases. 
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We begin with a more scrupulous analysis of the nature of 
the superconducting state. Since the character of the phase is 
determined by the dispersing fermion band, we will focus on 
the mean-field Hamiltonian for the down spins. If we com- 
bine the down-spin fermions on the two sublattices into the 
following spinor, 



fqU 
fq2i 



(14) 



then the Hamiltonian has the general form: 



down — 



e ix) u 



Jz (~ _ J_ 
Jz 



(*t g ) 



h.c. 



where we have taken J x = J v = J, and 



F (x) = 

Ay) . 



J 



\/3q y 



Qx . V ^Hv 

cos — sin - 

2 2 

_(l + 2cos _ cos _ 



J 



represent the kinetic energy for fermions hopping on the hon- 
eycomb lattice. The third line corresponds to an in-plane 
'magnetic field' in pseudospin space due to the enhanced hop- 
ping along the z links. This term shifts the positions of the 
Majorana cones, but is otherwise unremarkable. 

The second line is a superconducting pairing term along the 
x- and ?/-links. Both 
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J \Z3q y q x 
— cos — — — cos — 
8 2 2 



J . \f3q y 



•iy Qx 

sin — — - cos — 



(17) 



are non-vanishing in the mean-field state. The superscipts (s) 
and (i) refer to the fact that these are pseudospin-singlet and 
pseudospin-triplet superconducting order parameters. 

If we linearize about the nodes (we work at the isotropic 
point, J = J z , for simplicity), then the Hamiltonian for down- 
spins takes the form: 



H, 



down 



32 
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Vxj-^y 
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16 
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*>y(*- P ) T + h.c. 
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I [PyVx ~ PxHy] (* L p r + h 'C- (18) 



Here, p is the distance from the node (4tt/3, 0). This Hamil- 
tonian has four eigenvalues, the two non-dispersing ones 
± Jz/8, and the two dispersing ones in Eq. [13] 

It is helpful to isolate the dispersing band (unlike the Hamil- 
which contains both the dispersing and non- 
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toman 

dispersing down-spin bands). To this end, we form the Dirac 
fermion 



'/</ 



i7r/4 



(Cgl - ic q2 ) 



(19) 



The mean-field Hamiltonian for rj q is (up to a constant): 

^ = 9 J2 (Wll* + + h - C ) (20) 



where 




Qx v3 
2J cos — cos —^Qy 



1 Qx . 

A q = -Jcos — sm—q y 



(21) 



(22) 



To understand this Hamiltonian better, it is useful to mo- 
mentarily imagine that A q = and focus on e q . The Hamil- 
tonian now describes spinless fermions on the honeycomb lat- 

(15) 

tice with dispersion e q . First consider J z > 2 J. We see that 
there is no Fermi surface: e q is never equal to zero. Consider 
the minimum energy excitation, which occurs at q = (0, 
and has energy J z — 2 J. Near the minimum the band is ap- 
proximately quadratic. There are no excitations near zero en- 
ergy because the effective 'Fermi energy' lies below the bot- 
(16) torn of the band. Superconductivity does not change this pic- 
ture very much, other than to break U(l) symmetry (which is 
very important when we go beyond mean-field). When super- 
conductivity is turned back on, there are no nodes or nodal 
excitations because there is no Fermi surface. 

For J. < 2 J, there is Fermi surface which surrounds the 



point (0, ^t=). Strictly speaking, for the usual Brillouin zone 
this point sits on its boundary, so half the Fermi surface en- 
circles (0, ^7=) while the other half encircles the equivalent 

point (0, — ^) which differs by a reciprocal lattice vector. Of 
course, we could take a different unit cell for the reciprocal 
lattice which only includes one of these two points; then the 
Fermi surface will surround this point. We now restore the su- 
perconducting gap A q . This opens a gap on the Fermi surface, 
except at the points on the Fermi surface which intersect the 
nodal line q y — ^= ) . (The nodal line q y = does not inter- 
sect the Fermi surface, except for the point (47r/3, 0), which is 
equivalent to (27r/3, 2tt/ y(3)) under translation by a recip- 
rocal lattice vector.) For 2 — J z /J <C 1 small, the Fermi sur- 
face is approximately circular. Let us expand momenta about 
(0, so that {q x ,q y ) a (0, ^) + (2p x , 2 Py /V3). Then 

e p w J(p x +P%) ~ M> wnere the 'Fermi energy' [i is given by 
/i = 2 J — J z , and A p = J p y . Thus, the Hamiltonian in the B 
phase looks like that of a p y superconductor, which has nodes 
at Py = 0. As J z is decreased and the system moves towards 
the isotropic point, the nodes move towards the corners of the 
Brillouin zone, eventually reaching the graphene spectrum at 
the isotropic point. 



HI. MEAN-FIELD PHASE DIAGRAM IN THE ABSENCE 
OF TIME-REVERSAL SYMMETRY-BREAKING 
PERTURBATIONS 

We will now apply the mean-field description outlined in 
the previous section to understanding the phase diagram of ([T]) 
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in terms of its fermionic band structure and superconducting 
gap. As we shall see, the principle advantage of the spinful 
mean-field decoupling is that it allows us to better understand 
the system's behavior away from the exactly solvable point - 
both in terms of proximate phases, and the fate of physical 
quantities such as the spin-spin correlation functions as we 
perturb the Hamiltonian ([l}. At the end of this section we also 
describe at mean-field level the nature of the phase transition 
separating the gapped A phase and gapless B phase. 



A. The A Phase 

We begin by studying the A phase, for which J z > 2J 
and the band structure ( p"3j ) is fully gapped. In this phase, su- 
perconductivity, which couples fermions along the x- and y- 
links, competes with dimerization along the z-links, as is ev- 
ident from the 2-band Hamiltonian ( p"5j ). In the A phase the 
dimerization term dominates, leading to a fully gapped band 
structure, in the extreme limit J = 0, J z 7^ 0, dimerization 
leads to a gap, even in the absence of superconductivity. (In- 
deed, many fruitful explorations of the A phase treat it as an 
effective theory of such interacting dimers^Kil). 

As seen at the end of the previous section, we may view 
the A phase as a spin-polarized p-wave superconductor with 
chemical potential which lies below the conduction band. One 
amusing consequence of this is that the topological order of 
this phase is, as explained in Ref. Q3] that of a Z2 gauge the- 
ory. Its topological nature stems from the fact that, in the con- 
densed phase, the only remnant of the interactions between 
gauge fields and matter is a 'statistical' interaction due to the 
Berry's phase of ir accrued by a charge if it encircles a vortex 
of flux This provides an alternative perspective on the 

well-documented fact 1 11 that the A phase is smoothly con- 
nected to the so-called Toric code^ - a model of Ising spins 
which realizes a topological Z 2 gauge theory with matter. In 
particular, this highlights that the topological order of the A 
phase is not restricted to the set of exactly solvable Hamilto- 
nians described by ([Tji, but is that of a garden-variety s-wave 
superconductor. 

If we only cared about the single-particle gap, then we 
could close the superconducting gap entirely without closing 
the total fermion gap. However, the gauge symmetry of the 
problem would not be broken down to Z2 in this case, so there 
would be gapless gauge field fluctuations about the mean-field 
solution. (In the dimerized limit J x = J y = 0, though the 
£7(1) gauge symmetry is unbroken these gapless modes are 
absent since the gauge field cannot propagate). 

Because the A phase is fully gapped, it is stable to weak per- 
turbations away from the solvable point discussed here. For 
instance, we could add a weak magnetic field and/ or Heisen- 
berg interaction without changing the qualitative features of 
this phase. Since the system is fully-gapped, perturbation the- 
ory can be used, and the effect will be small, so long as the 
perturbation is weak. This is in contrast to the B phase which, 
as we will see, is unstable in the face of appropriately chosen 
perturbations. 



B. The Nodal B Phase 

We now briefly describe the B phase, for which J z < 2J. 
Now ( |20l i is the band structure of a p-wave superconductor 
whose nodes intersect the Fermi surface at two distinct points 
in the Brillouin zone. 

To simplify the algebra, we will consider the symmetric 
point J x — J y — J z = J. The energies of the dispers- 
ing Majorana bands are then exactly those of free fermions 
in a honeycomb lattice. The spectrum is gapless at the 
points q = (±^, (and at the equivalent points (±4f , 0), 



= ), which differ from the first two by reciprocal lat- 



v 3 ' ^> 

tice vectors). These nodes account for two distinct cones 
in the energy spectrum, as in graphene. However, unlike in 



graphene, the band structure ( 13 1 is that of a pair of bands of 



dispersing Majorana fermions. In the vicinity of these nodal 
points, it is useful to rewrite the Hamiltonian ( |20| in terms of 
the spinor 



Xq 



(23) 



where q is restricted to lie in half of the Brillouin zone to avoid 
double-counting, e.g. over q x > 0. In terms of this spinor, the 
Hamiltonian can be written in the form: 



H = O 4 ^1 Tx + e l Tz \ Xq 



(24) 



9x>0,q„ 



In the vicinity of the nodes (at the isotropic point J z = J), we 
can expand q = (^-,0) + (p x ,p y ) and write 



X P 



-(^,o)-#, 



(25) 



and p now ranges unrestricted over small p (e.g. over \p\ < 
A, for some cutoff A), i.e. near the nodes. Expanding e = 

^^-p x , we can write: 



V3J n a — V^ J i 



16 



H 



E4 



~32 



-Py 1 



~32 



~PxT z 



X P 



d 2 x x* U9, 



yTy + id X T z ] X 



(26) 



with v = J. Thus, these two Majorana fermions combine 
to form a single Dirac fermion. This Dirac cone is formed by 
combining the two nodes of a p y superconductor. This sin- 
gle Dirac cone does not violate the usual fermion doubling 
arguments since the gauge symmetry is broken. We will see 
presently, however, that it is central to the non-Abelian statis- 
tics of the gapped B* phase. 

We now consider some of the correlation functions of the B 
phase. Since there are gapless excitations, the energy-density 
will certainly have power-law correlations. How about the 
spin-spin correlation function? At the soluble point, this is 
short-ranged. Consider, for instance, the S z — S z correlation. 
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In terms of the slave fermions S. 
up and down-spins decouple, 



(4/*T-/i/4)/2. Since 



The first term vanishes since it only involves b x and b v , 
and these create/annihilate fermions in the up-spin flat bands. 
Here, b x and b v are defined in terms of /j,, /J according to 
Eq. 10 (It is important to remember that, although they play 
the same role in our analysis as the operators with the same 
labels in Ref. 1 they are not identical, in spite of the obvious 
similarity.) Thus, we are left with 



= (l + iibtc^ + iibfc) 
= 



<6f 



b~c 



,»/16 
(28) 



At the mean-field level, this is a free fermion problem, so 
we can evaluate these correlation functions. The Hamiltonian 
does not mix b z with c, so (ibfci) — and (bfcitfcj} — 
(bfbj) (cjCi). Since b z creates a fermion in a flat, non- 
dispersing band, (bfbj) = unless i and j are the same or 
neighboring sites. 

One of the appealing features of the formalism we use is 
that correlation functions in the presence of small perturba- 
tions to the Hamiltonian ([T} can be calculated with relative 
ease. For instance, suppose we consider a weak magnetic field 
in the z-direction, as in Ref. [T7] This adds a perturbation to 
the Hamiltonian: 



Hp eI t — o^zy^(/jf/»t _ filfil) 



(29) 



For small h z , this perturbation does not spoil the basic struc- 
ture of the spectrum: there are still three gapped bands and 
one gapless one. The up-spin gapped band will still be non- 
dispersing and will be at the same energy, but the correspond- 
ing eigenoperators will mix b x and b v (unlike the eigenop- 
erators (Hi in the unperturbed Hamiltonian). The down- 
spin gapped band will now disperse weakly, but will remain 
gapped. However, the eigenoperators for the down-spin bands 
will now mix b z and c. Thus, when we compute the (SfS*} 
correlation function, b z will have a small amplitude, propor- 
tional to h z for small h z , to create a dispersing fermion. Thus, 
this correlation function will have power-law falloff. 

To see this more precisely, we add the magnetic field term 
to the down-spin Hamiltonian: 



H. 



down 



JV3 

~32 



Jy/3 

~32 



J 



1G 
-h.c 



- ±*l[i y (*l p ) T + h.c. - lh z tf p % (30) 

When we diagonalize this Hamiltonian, we find a new set of 
eigenoperators a z ±, &o±. The eigenoperator a z + creates a 



fermion in a weakly-dispersing gapped band and has short- 
ranged correlation functions. The eigenoperator <So+ creates 
a fermion in a gapless band and has power-law correlation 
functions. For small h z (and, for simplicity, small momen- 
tum k), we can express the fermions a z ± = (ib z l ±b z q2 ) 1% 
a a ± = {ie l9q c q ^i ± c 9i2 ) /2, in terms of these new eigenop- 
erators as: 



at z ± 
a Q± 

Thus, we now have: 



a z ± ± y"o± 



(31) 



(b z b z ) = - ((ot z +,i + a z ^ il )(a z+J + a e -j)) 

= -({.a>z+,i + & z -,i + h z (a 0+ j - a _ il )/2) x 
{otz+,j + oig-j + h z (a Q+d - a _,j)/2 



+ ~A ((«o+,iao+,j> + (ao-.tao-j)) 



(32) 



Here, we have assumed, for the sake of concreteness and sim- 
plicity, that i and j are on the 1 sublattice. From the Hamilto- 
nian (30 1, we have for large separation |x — y| and to zeroeth 



order in h z : 



(ao+, x ao4 



) + (a _ )X ao-,y) = 
dw d 2 k ^(k y +ik x /2) e » k '( x -y) 
2tt (2tt) 2 



(33) 



Therefore, at long distances, 



(ao±, x ao±, y / 



(34) 



Combined with the (cjCj), which has the same-power-law, this 
gives an (S Z S Z ) correlation function which falls off as 1/r 4 
in the presence of a small magnetic field, in agreement with 
the results of Ref. [T7] 

In the face of perturbations that are not quadratic in the 
fermions, such explicit calculations are more difficult in gen- 
eral. However, as is frequently the case in spin-liquid models 4 , 
the structure of the Fermi surface (here a pair of Dirac cones) 
is protected by symmetries of the mean-field state. Thus small 
perturbations which do not break any symmetries of the prob- 
lem cannot open a gap in the spectrum. 



C. Transition between A and B Phases 

As we move within the gapless B phase, from the isotropic 
point J x — Jy — J z towards the boundary to the A phase, the 
two nodal points move together and, at the phase transition 
point, merge. The nodes then annihilate as the phase boundary 
is crossed. In this section, we focus on the transition point. 
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As discussed in Section [TTJ the dispersing spin-down band 
can be rewritten as a model of spinless fermions with p y su- 
perconducting order, as in Eq. 20 At the boundary be- 



tween A and B phases, the Fermi surface has shrunk to a 
point because the effective chemical potential is precisely at 
the bottom of the band. When the effective chemical poten- 
tial is at the bottom of the band, the spectrum is quadratic 
in the absence of superconductivity. Superconductivity with 
p y pairing symmetry leaves the spectrum gapless but makes 
the spectrum linear in one direction. We now examine this 
in more detail. Expanding about the bottom of the band 
(Qx,Qy) ~ (0, + (2p x ,2p y /V3), we can write the Hamil- 
tonian d20l in the form: 



H = 



2 ^ 

p 



g P 2 vtVp ~ jPy {vW-p - VpV-p) 



T 

X—p 



Px>0,p y 



where 



rPv 1 ~ nP 2iT v 



Xp 



If we go to a Majorana basis, 



this can be re-written: 



V P + V-p 



V2\(v P ~ V-p)/iJ ' 



(35) 



(36) 



(37) 



H 



1 E 

2 ^ 

Px>0,p B 



T 
<P-p 



J J 



P T y 



•rv 



d 2 x <f 



-id T - -8 2 t 
4 LU v 'z 8 



'■P 



(38) 



Therefore, the low-energy theory can be called a single 
gapless Majorana fermion, albeit an anisotropic and non- 
relativistic one. 



IV. BEYOND MEAN FIELD THEORY 

Thus far, we have found a consistent mean-field solution 
of ([T]i using the fermionization <J3J which reproduces exactly 
the Majorana fermion band structure and phase diagram of the 
exact solution proposed by Ref. 1 , We next ask what can be 
said about its fate upon including fluctuations of the various 
bosonic fields. The answer is not obvious since, unlike the de- 
coupling used by Ref. 1 ] the product bfbf +1 on each link does 
not commute with the full unprojected fermion Hamiltonian 
(although it does commute with the quadratic Hamiltonian 
Hmf)- Here we first establish that these fluctuations do not 
alter the results of the previous sections. Second, we demon- 
strate that at long wavelengths these bosonic modes lead to 
precisely the Z2 gauge theory of Ref. 1 Together, these facts 
cement the equivalence between the fermionization (|2| and 
Kitaev's exact solution. 



The underlying reason for this stability is that the unpro- 
jected mean-field wave functions we obtain can be mapped 
via Eq. (10 1 onto unprojected wave-functions in the Majo- 
rana fermionization of Ref. [T] Enforcing the SU(2) gauge 
constraints to reduce the model back to the physical Hilbert 
space amounts to two things: first, it eliminates the disten- 
tion between different possible mappings between f a . /J and 
b x .y,z,c. Second, it imposes a condition which is equivalent 
to the Z2 constraint required for the fermionization of Ref. 
1 Thus when expressed in the Majorana basis given by ( [T0| , 
the effect of this projection will be to apply the projector rel- 
evant to Kitaev's Majorana fermionization. In this way, both 
fermionizations lead to the same wave functions after projec- 
tion. 



A. Symmetries and robustness of the mean-field solution 

First, we will show that for the solvable Hamiltonian ([T}, 
the model's unusually large number of symmetries protect the 
exact fermionic band structure. The mean-field solution is 
thus exact, in that it correctly describes all correlators of the 
physical spin degrees of freedom, in spite of the apparent vio- 
lence done to the wave-function by Gutzwiller projection. 

We begin by listing the symmetries which are relevant to 
this discussion. The Hamiltonian ([T]i has the following dis- 
crete symmetries 



-tx,y,z 



1 <? : 



x,y,z 



(39) 



where the sign s XtV , z — ±1 can be chosen independently for 
x, y, and z spin operators. In the fermionic description, this 
leads to two discrete symmetries preserved by the mean-field 
Hamiltonian: 



C 

s 



fqai f-qo-i 



-q2a 



(40) 



Here, the 'charge conjugation' symmetry C is unitary (it is is 
simply i' qiaa -> ir x ) ab ip qiab ) while the 'sublattice' sym- 
metry S is an anti-unitary symmetry. Thus, in the mean- 
takes A,j , tij — > 

Quadratic Hamiltonians invariant under 



field Hamiltonian, C takes A^- , tij — > Ay , tij while S takes 



A • ■ t ■ ■ — v A* t* 



C have eigenstates which are diagonal in the Majorana basis 
( 10 1. These symmetries impose an important restriction on tij 



and Aij. C is preserved as long as ty, Ay are purely imag- 
inary. S is preserved so long as there are no terms directly 
coupling fermions on the same sublattice. 

Time-reversal symmetry is also respected by the model and 
its mean-field solution: 



9,t 



J—QiX 7 fq,-l ^ f—q. 



t 



(41) 



Single-spin terms (i.e. a magnetic field) and three-spin in- 
teractions break this symmetry. However, not all T-breaking 
perturbations will open a gap in the B phase: only those per- 
turbations which break S will open a gap in the spectrum, as 
we will see below. For example, the magnetic field discussed 
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in Sect. MB breaks T and C, but not S. As shown explic- 
itly above, this does not gap the B phase and indeed results in 
power-law spin-spin correlations. 

The relation between these symmetries is: 



SG X C 



(42) 



where the symmetry G x is given by: 



(43) 



which are a discrete subset of the off-diagonal SU(2) rota- 
tions interchanging up and down spins. In the mean-field so- 
lution these are no longer local symmetries. However, they 
remain global symmetries of the theory, whose effect is to 
rotate between different possible mappings between the four 
Majorana fermions (c, b x ' v ' z ), and the four self-adjoint combi- 
nations fj a + fi a , i (flu — /jo- J of the spinful fermions. Thus 

T is a projective symmetry - a symmetry that maps the sys- 
tem to a different but gauge equivalent saddle point. Such 
projective symmetries are important to classifying the phases 
of spin-liquid systems^. 




FIG. 1: [Color Online] The product of spin operators conserved 
separately on each plaquette by the Kitaev Hamiltonian {T]l. The 
Hamiltonian (TTj distinguishes between three types of links on the 
honeycomb lattice, which we call x-, y-, and z- links ( color-coded 
red, green, and blue respectively here). On x-links the spin-spin in- 

(x) (x) 

teraction term is Sj , and similarly for y- and z- links. The 
product of spin operators shown here - a product around a plaquette 
of the spin variable associated with the 'external' edge at each vertex 
- commutes with the spin Hamiltonian. 



Besides these more generic discrete symmetries, Eq. ([T]i 
represents a somewhat special point in a more extended space 
of similar spin Hamiltonians: there is a product of spin opera- 
tors on each plaquette which commutes with H. This is: 



2 6 



(44) 



where e(i) — z for a vertex which sits between x and y links 
on the plaquette, y for a vertex which sits between x and z 
links on a plaquette, and x for a vertex which sits between y 
and z links on a plaquette (see Fig. [IJ. In the ground state, the 
value of this operator is positive on each plaquette 1 . 
In terms of the fermionic operators, V can be written as 



(45) 



where Pq denotes Gutzwiller projection onto singly occupied 
states, a = x,y, z on x,y, and z-links, respectively, and bf 
are the Majorana fermions defined in Eq. (10 1. (Since the 



quantity in parentheses is not SU(2) gauge invariant, the pro- 
jection operator is necessary in this case). In the mean-field 
state, each species of Majorana fermion is localized on the ap- 
propriate links, with (bfbf +1 ) mf = 1/2- Terms annihilated 

by P do not contribute, since (f^f-^MF = (MM)mf 
0. Hence we find that the mean-field value 



V f = (blbl)(blbl)(blbl)(blbt)(blbl)(blbl) 



1 



(46) 



is precisely that of the exact solution. 

We now show that, combined with the discrete symmetries 
mentioned above, conservation of Vf prevents fluctuations 
about mean-field from altering the fermionic band structure 
in any way. We will first establish that the symmetries for- 
bid any terms other than those in Eq. (ftp} from contributing 



toVf. If there can be no further contributions to Vf induced 
by fluctuations, however, then also no spectral weight can be 
transferred from the equal-time correlation functions of the 
b a , as otherwise we would not arrive at the correct value for 
V . This means that all further-neighbor correlators must van- 
ish exactly. 

By Wick's theorem, we need only consider the possibility 
of other pairings of the fermionic operators which give a non- 
zero contribution to Vf. The only possibility allowed by C 
and PS is to give a non-vanishing expectation value to terms 
of the form (bf etc. Thus we consider: 



{blbt)(blbl)(blbl)(blbt)(b\bl)(blbl 



(47) 



However, the interacting Hamiltonian for the up spins decou- 
ples exactly into separate Hamiltonians for each chain of x — y 
links in the lattice. In particular, the full Hamiltonian con- 
tains no interaction term coupling bf and b\ , as they lie on 
different chains. Hence interactions cannot shift (bfbf ) from 
its mean-field value of 0. As (47i is the only extra contribu- 



tion to Vf not explicitly forbidden by symmetry, we conclude 
that Eq. ( |46| must remain valid in the full solution, and that 
consequently no fermion bilinears can be shifted from their 
mean-field values. 



B. Gauge theory of fluctuations about mean field 

Thus far, we have shown how to reproduce Kitaev's mean- 
field portrait of the exact spin-liquid ground state using the 



9 



fermionization Q, and argued that including fluctuations 
about mean-field will not change the fermionic band structure. 
Hence we have obtained an alternative mean-field description 
of the ground state of ([TJ which reproduces faithfully the spin 
correlators of the exact ground state. 

Though the mean-field solutions describe identical physics, 
however, the fermionization (|2j differs quite dramatically 
from that of Ref. [TJ in the nature of the bosonic variables, 
and consequently, the theory of fluctuations about mean-field. 
After Hubbard-Stratonovich transforming the 4-fermion inter- 
actions, we obtain bosonic fields which condense to give both 
the hopping and superconducting order parameters, as well as 
the SU(2) gauge fields associated with the constraint pj. One 
might therefore wonder why these do not lead to significantly 
different physical theories after fluctuations about mean-field 
have been accounted for. Here we address this question, al- 
lowing us to posit that (|6j describes a gapped spin-liquid phase 
which exists even away from the exactly solvable limit of the 
Hamiltonian dT). 

The bosonic fluctuations about mean-field can be separated 
into the following degrees of freedom. There are three scalar 
fields describing fluctuations in the amplitudes of the various 
kinetic and superconducting terms. All of these are massive, 
and as we shall see two of them can be interpreted as Higgs 
fields for the broken SU(2) symmetry. In addition, there 
are three independent fields associated with phase fluctuations 
of the various link variables. These can be identified as an 
SU(2) gauge field (describing phase fluctuations of the spin- 
symmetric hopping term) and two 'Goldstone bosons' asso- 
ciated with the phases of the order parameters breaking the 
SU(2) symmetry. We will briefly discuss each type in turn; a 
more detailed analysis is presented in Appendix |B 3| 

We begin with the scalar fields describing fluctuations in 
the amplitude of the various bosonic order parameters that fix 
the mean-field fermionic band structure. The general form of 
the Hubbard-Stratonovich action ensures that all of the scalar 
fields are massive, with energy gaps of order j at the isotropic 
coupling point. Because of this mass gap, fluctuations in the 
amplitudes of the mean-field parameters are not generally ex- 
pected to have an important effect on the fermions. The no- 
table exception to this^ is in cases when they destabilize the 
spin liquid saddle point in favor of a 'dimerized' state with 
spins hopping predominantly along a subset of links in the 



lattice. As we discuss in Sect. Ill A an analogue of the dimer- 
ized phase does occur for anisotropic J x ,y,z'< m general we 
may therefore conjecture that away from the solvable point 
this phase boundary may be shifted, but that fluctuations of 
the mean-field hopping and superconducting amplitudes will 
not qualitatively alter the phase diagram. 

Next, we consider the impact of phase fluctuations de- 
scribed by the SU(2) gauge theory. Naively, the gauge theory 
is strongly-fluctuating, since there is no small parameter in the 
problem. However, the ground state of ([TJ is a Higgs phase, 
so that the gauge field is massive. (Importantly, this explains 
why the gauge theory is not confined). 

To see that the model ([TJ is in a Higgs phase, we view the 
mean-field solution ([6]) as a condensate of two independent 
order parameters in the adjoint representation of SU(2). As 



explained in detail in Appendix B 3 the combination of su- 
perconducting and spin-antisymmetric hopping terms break 
the SU(2) gauge symmetry. This leaves only the residual Z 2 
gauge symmetry group one normally finds in a superconduc- 
tor: 



fieri fia ^ fieri f i a ^ij.o , ^-ij,cr 



ft 



(48) 



comprising the residual Z 2 symmetry of the U{1) subgroup 
broken by superconductivity. As a result of the Anderson- 
Higgs phenomenon, the dynamical fluctuations in the gauge 
field are suppressed at long wavelengths, so that gauge 
field fluctuations are not expected to substantially alter the 
fermionic band structure. (Here the gauge field results from 
the constraints of the purely 2 dimensional system, and conse- 
quently is fully gapped unlike the electromagnetic gauge field 
in thin-film superconductors.) However, the gauge field makes 
itself felt in the interesting topological structure of the spin- 
liquid phase. 

An alternative route for a gauge field to acquire a mass is 
through the generation of a Chern-Simons term. We will re- 
turn to this possibility when we consider perturbations break- 
ing T in Sect, [vj where we shall see that it plays an important 
role in the topological nature of the theory. 

In summary, we can understand the exact ground state of 
([TJ - a phase whose propagating degrees of freedom consist 
of Majorana fermions coupled to a Z 2 gauge field - as a rather 
special incarnation of the Z 2 spin liquid: a spin-polarized p- 
wave superconductor. In this description, we arrive at Majo- 
rana fermions not by expressing the spins directly in a Majo- 
rana basis, but rather by starting with Dirac fermions coupled 
to an SU(2) gauge field and choosing a mean-field solution 
which breaks the gauge symmetry. The Z 2 flux is thus the 
superconducting vortex, while the Z 2 charge carried by the 
Majorana fermions reflects the fact that the superconducting 
state conserves charge modulo 2. 



T-BREAKING PERTURBATIONS: THE GAPPED B* 
PHASE 



In the previous section, we showed that one way to open 
a gap in the B phase - by merging the two nodes - can be 
understood as a transition between a nodal and nodeless su- 
perconductor. This drives the system into the A phase. There 
is, however, a second way to open a gap: we may add an- 
other pairing term to the effective Hamiltonian ( p"5j ), which 
will fully gap the spectrum provided that the corresponding 
gap does not vanish at the Dirac points. Here we focus on this 
latter gapped phase, and discuss its topological properties. 

As noted in Sect. |IV A| this second gapped phase neces- 
sarily breaks one of the two discrete symmetries of the mean- 
field solution - and hence the physical time-reversal symmetry 
of the spin model - since we must include couplings between 
sites on the same sublattice. Here we will focus on the case of 
broken S, as this can be realized by adding a 3-spin interaction 
which commutes with the Hamiltonian ([TJ. 
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A. Mean-field theory with T -breaking terms 

In terms of the original spin degrees of freedom, the T 
breaking term we must add to enter the B* phase is: 



.r_ 

2 



£ 



s 



; qz qV 



qz qV qx 
°i °j °k 



£ 



QZ qx qV 
°i °i J fc 



(49) 

(see Figure [2]). It is easy to see that this commutes with the 



plaquette product of spins (44, , and hence preserves the Z2 



vorticity on each plaquette. Hence it also commutes with the 
full Hamiltonian- though not individually with the spin bilin- 
ears on each edge. 




FIG. 2: [Color Online] The 3-spin interaction which breaks T but 
commutes with the conserved 6-spin product on each plaquette (see 
Fig. [T|r 19 . The operator is constructed by taking the product of spin 
operators on three adjacent vertices, where the direction of the central 
spin is that associated with the 'external' edge at the vertex, while the 
two external spins match the edges joining their associated vertices 
to the central vertex. In the figure the central vertex is at the upper 
left, and the operator is the product of a z at the lower left vertex, a y 
at the upper left vertex, and a x at the top vertex. 



Expressing the spins in terms of Dirac fermions yields a 6- 
fermion interaction. Though we cannot perform the analogue 
of an exact Hubbard-Stratonovich transformation for the re- 
sulting action, which contains both 4 and 6 fermion terms, at 
small J' it is possible to evaluate its effect on the mean-field 
solution in a controlled way (see Appendix [C]). We find that 
(consistent with the treatment of Ref. [T|) the effect of such 
a term is to induce second neighbor hopping and supercon- 
ducting terms, without altering the rest of the band structure 
(except for an overall rescaling of the bandwidth). 

We therefore begin by studying the resulting mean-field 
Hamiltonian. The 3-spin interaction introduces the following 
quadratic fermion terms for the down spin band: 



H 



(i) 

MF 



J' 



= —(-sinq x +sinq- W - sing - Z 2 ) [-^\p. z ^ q 

+m^_ q ) T +h. c ] . (so) 



where "J was defined in Eq. 14 As shown in Appendix [C] 
the perturbation ( |49| ) does not alter the mean-field Hamilto- 
nian of the up spins, which therefore maintain their flat band 
structure and remain localized on x- and y-links. In addition, 
the new couplings do not disrupt the pair of flat spin down 
bands. Thus the basic structure of the initial mean-field solu- 



tion is preserved, and the only effect of the interaction ( 49 1 at 
mean-field is to alter the structure of the dispersing spin-down 
band. 

The new effective mean-field Hamiltonian for the spin- 
down fermions therefore has the form: 



7 - <-i J \ fl/t 

J, 



*J/Xy*, (51) 



with e q x ' v \ Aq S ' p - ) given in Eq.s ( 16|l7 1, and 



J' ( . . „ . Qx V3q y 
— I - sin q x + 2 sin — cos 



(52) 



In the vicinity of the Dirac cone, for J x ,y,z = J, this gives: 



down 



Q 



V3 T v/3 J 

J'q 2 - ' -J' \ H 



SVS T/ _a 3^3 A 
-j 1 /i 



\ 64 16 

+ l^^P [PyVx-PxVy] (*l p ) T +h.C. 

-^^(*i P ) T + h.c. 
+ (3/8V3J' 2 g 2 + 3^3J' /2)^\{^_ q ) T + h.c. (53) 



which we can view as a mixed s- and chiral p-wave supercon- 
ductor. This term opens a gap at the Dirac cone, so that the 
system is now fully gapped. We discuss the consequences in 
the next subsection. 



B. Topological features of the gapped B phase 

Thus far, we have established that adding the spin interac- 
tion (49 1 has the effect, at mean-field, of breaking S and open- 
ing a gap in the spectrum of the dispersing Majorana mode 
(c), whilst leaving the band structure of the localized Majo- 
rana modes (b x ' v ' z ) unchanged. We will now see how this 
perturbation leads to a topological phase with 0-energy Ma- 
jorana fermions bound to vortices, exactly as in the spinless 
p + ip superconductor of Read and GreerPn 

The simplest way to identify the nature of the B* phase is 
to consider the Hamiltonian d20|, where the B phase is a p y 
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superconductor. The perturbation modifies the Hamiltonian 
according to: 



A„ — > A„ — i— sin 



.J' . q x ( \/3 



cos q v — cos ■ 



■isgn( gx )^(l + §) yjl-{^f (54) 



In the second line, we have approximated A 9 by its value in 
the vicinity of the nodes. From this expression, we see that 
this is an ip x superconducting gap which opens up a gap at 
the nodes. 

As noted previously, in the nodal B phase, the 'chemical 
potential' p, = 2 J — J z lies in the band. Thus, when the gap 
is opened, the system goes into the 'weak-coupling' p + ip 
superconducting phase. As q ranges over the Brillouin zone, 
the vector (ReA g ,ImA g , e q )/(e q + lAgl 2 ) 1 / 2 wraps around 
the sphere. The corresponding winding number cannot be 
changed without closing the gap, i.e. without going through a 
phase transition. 

Conversely, when the 3-spin interaction is included in the 
A phase, the chemical potential lies below the band. For 



sufficiently small J', (ReA g , ImA 9 , e,)/(e 2 + lA^ 2 ) 1 / 2 re- 
mains in the northern hemisphere, and thus has winding num- 
ber zero. Thus, this is the strong-pairing phase of the chiral 
p-wave superconductor. In other words, including a weak Si- 
breaking perturbation in the A phase leaves the system in the 
A phase. 

Once we have identified the B* phase with the weak-pairing 
phase of the chiral p-wave superconductor, we are faced with 
the following riddle: in its usual incarnation, the supercon- 
ducting coherence length is assumed to be much larger than 
the lattice scale, so that vortices are well-modeled by a con- 
tinuum theory. In particular, the vortex will have a core which 
is in the normal state. The argument put forth by Read and 
Green 2 ^ to show that in the weak-pairing phase a 0-energy 
Majorana fermion is bound to the vortex core relies on the 
existence of a 'domain wall' between the vortex core and the 
superconductor in an essential way. Since phase B* is known 
to have the same topological order as the chiral p-wave super- 
conductor, in which the existence of Ising anyons is due to the 
fact that these 0-energy Majorana fermions are bound to the 
vortex cores, we expect a similar phenomenon. In the lattice 
model at hand, however, a vortex exists on a single plaque- 
tte, and there is no vortex core. How, then, do the Majorana 
fermions become bound to these vortices? 

One answer to this question comes from studying the long- 
wavelength gauge theory. First, we observe that the key effect 
of the T-breaking 3-spin interaction is that it induces a mass 
term -, = — to^^t Q s = ^^J' at the two nodes in 

the Brillouin zone. As discussed perviously, the low-energy 
effective theory is that of a single species of massive Dirac 
fermion. If we integrate it out, then as shown explicitly in 
Appendix [P] the 1-loop effective action for the gauge fields 
is precisely what we would expect from a single Dirac cone, 
except that, since U(l) is broken down to Z 2 , a Higgs mass is 



also generated: 
£ (iloop) 1 2 



Airm 
to 1 



K A' i d u A > 



(55) 



In other words, we obtain the usual Higgs mass term, the field- 
strength tensor squared, and a Chern-Simons term with level \ 
(as usual from a single Dirac cone^j). The Higgs mass is pro- 
portional to the condensate fraction |$| 2 , and is crucial outside 
a vortex. However, in a vortex core, the condensate vanishes. 
We will assume that the Higgs mass can be neglected in the 
core. Thus, in a vortex core, we have 



SC 
8A^ 



TO 87T 



(56) 



where is the fermion current, and we have used d v F llv 
0. Taking p = 0, to > 0, we obtain the constraint: 



1 

47T 



B n = pa 



where p = Jq. In the case at hand, we have 



Pq 



\jk,ifk-q,i + f-k,if-k+q,: 



i=l,2 k 



(57) 



(58) 



(Here k is technically restricted to momenta near the Dirac 
cone; more generally, we sum over only half the Brillouin 
zone.) The rather counter-intuitive fact that holes at the left 
Dirac cone carry the same charge as particles at the right Dirac 
cone results from the fact that the two cones have opposite 
chirality. 

The density p q is a sum of the density of particles at the 
right Dirac cone, and holes at the left Dirac cone. The creation 
operator associated with this density is the Majorana fermion 

c q = i f/tj — f-qij, which simultaneously creates a particle 

at q and a hole at — q. Hence Eq. (57i tells us that there is a 



Majorana fermion c bound to every half-flux quantum. These 
half-flux quanta are precisely the Z2 vortices of the supercon- 
ductor; hence we conclude that there is a Majorana fermion c 
bound to each Z 2 vortex. 



VI. SPIN-DENSITY WAVE STATES 

As described in the previous section, a T-breaking 3-spin 
term opens up a gap which can be written as follows in terms 
of the x fermions, 



X P 



Jl-Q/2-p; 



(59) 



At the isotropic point, J x = J y = J z , Q/2 = (^,0). The 
Hamiltonian in the B phase can be written in the form 



H 



[vPyTy + VP X T Z ] X P 



(60) 
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where v = ^p- at the isotropic point. The Dirac mass term 
generated by the three-spin interaction is of the form: 



Hum. = m^XpTyXp- 



(61) 



where m — 3J'/2. 

However, this is not the only possible term which can open 
a gap at the nodes of the B phase. The other possible term 
is (W is a coupling which we introduce to parametrize the 
strength of this term): 

■ffpair = W^2XpiTyX-p+h.C. 

P 

= 2W J2 r )U/2 + />$/2+p- + h - C - 
P 

= W E [ c Q/2~ p mCq/2 +p m + (1^2) 
p 

+ iC Q/2-p,l C Q/2+p.2 + (1 2 )] + h ' C - 
= - 4W/ 5Z[^/2-p,l^<3/2+p,l _ ^Qj2-pJ-Ql2-p,\ 



h.c. 



(62) 



Thus, such a mass term breaks translational symmetry. It 
includes terms which induce superconductivity at non-zero 
wavevector as well as terms which induce a spin-density wave 
at wavevector Q. We can imagine that a spin-spin interac- 
tion which is added to the Kitaev model as a perturbation 
will, upon decoupling, generate such a mass term. However, 
since the density-of-states at the nodes is zero, interactions 
will only generate such a term at O(l) coupling strength (not 
at infinitesimal coupling, as would the case for a Fermi sur- 
face instability). At O(l) coupling strength, there is no reason 
to focus on the the nodal regions, so many other instabilities 
could also occur. It is possible that, in a large-iV version of 
this model, such a translational-symmetry-breaking instabil- 
ity will occur at weak-coupling. 

Similar but distinct spin-density-wave states have recently 
been discussed in the context of a hybrid Kitaev-Heisenberg 
model in Refs. 121 1221 



VII. DISCUSSION 

In describing the spin-liquid ground states of the various 
phases of Kitaev's honeycomb model using the slave-fermion 
approach, we may learn several things about the nature of the 
phases of this model, their potential stability to perturbations 
away from the solvable point, and their precise relationship to 
other phases of matter which exhibit similar physics. 

First, the fermionic mean-field theory allows us to relate the 
various phases of the Kitaev model to the ground states of dif- 
ferent Bogoliubov-de-Gennes Hamiltonians. This can be done 
in two different ways: (1) in terms of the fermions /^j, intro- 
duced in Eq. [2] and (2) in terms of the fermions -q introduced 
in Eq. 19 The latter are formed from the propagating part of 



fl. Each way has its conceptual and technical advantages, as 
we have seen. 

The mean field phase diagram is summarized in Figure [3] 
which can be interpreted in terms of the 77 fermions as follows. 
The A phase, in which the nodes of the superconductor do not 
intersect the Fermi surface, is adiabatically connected to an 
s-wave superconductor. The B phase is a nodal p-wave super- 
conductor. The B* phase is the weak-pairing phase of a chiral 
p-wave superconductor, with the consequent Ising topologi- 
cal order. The A* phase is the corresponding strong-pairing 
chiral p-wave superconductor phase. As a result of the strong- 
pairing nature of this phase, the topological order is, in fact, 
again that of an s-wave superconductor. The reason for this 
is that, at the mean-field level (i.e. when treated as a free 
fermion problem) the A and A* phases can be adiabatically 
deformed into each other, so the line between them in Figure 
|3]is a crossover line. On the other hand, the other transitions in 
Figure [3] are genuine phase boundaries which are essentially 
the same as the corresponding transitions in the superconduc- 
tor. One important difference needs to be emphasized. In 
a two-dimensional superconductor with a three-dimensional 
electromagnetic field, there is a gapless plasmon. Thus, a 
thin superconducting film is not fully gapped, even though its 
fermionic spectrum is fully-gapped. However, in the Kitaev 
honeycomb lattice model, the gauge field is two-dimensional. 
Consequently, the plasmon is gapped and the system is fully- 
gapped. 



J' ± 



J/Jz 




p-wave super 
conductor 



FIG. 3: Schematic phase diagram of the Kitaev honeycomb model, 
and the corresponding superconducting phases. The phase is deter- 
mined by the ratio J/J z , and by whether the coefficient J' of the 
3-spin interaction is non-vanishing. 



Although the SU(2) mean-field theory described here is 
clearly more complicated than Kitaev's at the soluble point, it 
has the salient virtue that it is well-suited to perturbing away 
from the soluble point - partiaularly in the gapless B phase. 

It is interesting to consider the fate of the phase diagram 
shown in Figure [3] when the spin Hamiltonian is deformed 
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away from the exactly solvable point. The fully-gapped 
phases - the A and B* phases - will be robust against small 
perturbations by virtue of their energy gaps. So long as the 
gauge symmetry is broken to by the saddle point solution, 
the gauge field is gapped, and we do not expect fluctuations 
to lead to confinement. Therefore, the model still admits ef- 
fective spinon excitations, and the topological order of the 
spin liquid will be robust to gauge field fluctuations. Since 
the spinons are also gapped in these phases, they are stable 
against adding weak interactions between the fermions. The 
gapless B phase is a little trickier. Since the gauge field is 
fully gapped, we believe that the gauge field action is robust 
against small perturbations. The fermions, on the other hand, 
are gapless. However, since they have a single gapless Dirac 
point (rather than a Fermi surface), weak interactions between 
the fermions are irrelevant by power-counting. This is the rea- 
son that an 5£/(2)-invariant Heisenberg perturbation does not 
lead to a phase transition until the perturbation is sufficiently 
strong. Thus, we could say that the stability of the gapless B 
phase relies on phase space limitations. However, as we have 
seen, although the gapless B phase is stable against weak per- 
turbations, some features of the soluble point are not generic 
to this phase. For instance, a magnetic field will make the 
spin-spin correlation functions have a power law, rather than 
short-ranged, form. 

The fact that the bosonic fluctuations are all gapped does 
not, however, prevent the theory from acquiring a new lowest- 
energy saddle point if we deform far enough away from the 
solvable model. For instance, as we have discussed, the gap- 
less B phase can acquire a gap by an alternative method: the 
development of a spin-density wave, as discussed in Section 



VI 



Various perturbat ions of the Kitaev model, including a 
Heisenberg interactiorP ^jca n lead to such an instability. Fur- 
thermore, it is well known 3 18 that symmetric spin-liquid states 
are often prone to dimerization instabilities, in which the spins 
pair with neighbors in a valence-bond crystal which breaks a 
lattice symmetry. Away from the solvable limit, therefore, 
it is likely that the phase diagram will also include some 
such valence-bond crystal states. At the symmetric point, the 
model has a spin-orbit type 3-fold rotation symmetry (entail- 
ing a 3-fold lattice rotation about a vertex, coupled with a 
global spin rotation of the form ( |43"j which makes the saddle 
point perturbatively stable - though in principle lower-energy 
symmetry-breaking saddle points might exist. Away from the 
isotropic point J x = J y = J z , such states need not break 
any symmetries of the Hamiltonian, so that symmetry does 
not prevent the saddle point from flowing to such a valence- 
bond crystal upon including fluctuations of the amplitudes of 
the mean-field hopping and superconducting terms. 

The fact that the exact ground state of ([T]i can be correctly 
described in the slave-fermion mean-field approach used here 
is also interesting in its own right. As discussed above, since 
the mean-field state is a Higgs phase of the gauge field, the 
model is in a regime where the spin-liquid saddle point is most 
likely to be stable. Even in this case, however, examples of 
Hamiltonians where the exact ground state can be shown to 
be a spin liquid are rare. The Kitaev model is thus a poten- 
tial testing ground for the slave-fermion approach, since we 



may begin with a Hamiltonian for which it is demonstrably 
valid, and consider the fate of the ground state under various 
perturbations. In particular, on general ground^ we expect 
that for small perturbations which do not close the gap in the 
spectrum, the slave-fermion mean-field theory will continue 
to capture the topological order of the gapped phases. 

Another interesting prediction of the slave-fermion ap- 
proach is that near the solvable point, the Kitaev model be- 
comes a superconductor upon doping. Specifically, we imag- 
ine starting with a Mott insulator whose effective Hamiltonian 
at half-filling is given by (jlj. After doping away from half- 
filling we must account for the fermion hopping terms, lead- 
ing to a t — J model, with the spin Hamiltonian given by 
Following the prescription used to study the cuprates^, we 
may decompose the spin operators as in Eq. Q, and express 
the electron operator as 



flh 



with the constraint 



flh + flh + b\bt = 1 



(63) 



(64) 



It follows that, at temperatures below the Bose condensation 
temperature of the bosons, and at sufficiently low dopings that 
the mean-field solution described above is a good approxima- 
tion for the spinons (fi a ), the superconducting order parame- 
ter is: 



,phys 



( C /L C -fcer') — (fk+q,af-k-q,a')(b-q<jbqcr>) 
^k-a,cr'Ps (65) 



where p s is the bosonic superfluid density. Thus the momen- 
tum dependence of the physical superconducting order param- 
eter is set by that of the mean-field superconducting order pa- 
rameter A for the fermionic spinons /. For the Hamiltonian 
([TJ, this predicts spin-triplet superconductivity (with equal 
spin pairing), with a mixed singlet and triplet pseudospin or- 
der parameter. 

Finally, it is interesting to compare the mean-field ground 
state of the Kitaev model with existing proposals for gener- 
ating the B* phase's topological Majorana fermions in phys- 
ical materials. The mean-field Hamiltonian of the B phase 
is manifestly equivalent to a p + ip superconducting state of 
spin-polarized fermions^. It also has an interesting relation 
to the the effective Hamiltonian of Fu and Kana^for surface 
states of a topological insulator in the presence of induced s- 
wave superconductivity. In the absence of superconductivity, 
these surface states form a single Dirac fermion. This Dirac 
fermion is analogous to the Dirac fermion which we have in 
the gapless B phase. If a magnetic film is brought into con- 
tact with the topological insulator, and the magnetic moment 
is perpendicular to the interface, then the resulting term in the 
Hamiltonian is a Dirac mass term, which breaks time-reversal 
symmetry and opens a gap. This is analogous to the 3-spin 
term in the Kitaev model, which opens a gap and drives the 
system into the B* phase. Note that this term in the Kitaev 
model is not analogous to the term generated by an s-wave 
superconducting film on the surface of a topological insulator. 
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Instead, s-wave superconductivity on the surface of a topo- 
logical insulator is analogous to a term x^iT y x p +\\.c, which 
is a down-spin density wave at wavevector (87r/3, 0) at the 
symmetric point J x = J y = J z . 

In all cases, the essential ingredients for generating topo- 
logical Majorana fermions are a 2-band model in which the 
band structure is that of a massive Dirac fermion, and with in- 



duced superconductivity. As we described in Sect. V B above, 
the massive Dirac fermion in all of these models is implic- 
itly coupled to a gauge field, since it forms a superconducting 
state. The fermion mass therefore generates a Chern-Simons 
term in the effective gauge-field action, which has the effect 
of binding a half-quantum vortex to each charge, since there 
is only a single Dirac cone. The charge which is bound in 
the superconducting state is a Bogoliubov-de-Gennes quasi- 
particle, rather than a fermion - which, when the supercon- 
ducting order parameter has a p-wave component, binds a Ma- 
jorana fermion to the vortex. 
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Appendix A: Mapping between SU (2) and Majorana 
fermionizations 



Here we explain in more detail the correspondence between 
the fermionization (|2]i and the Majorana fermionization em- 
ployed by Kitaev 1 . We begin with the mean-field correspon- 
dence: 



b x 



b y - 
q/ 



/qit + f-lif 
Cqi = i (fqil — f-qil 



(Al) 



which gives a mapping between unprojected spinful fermions 
and unprojected Majorana fermions. This mapping is not 
unique, as each Majorana fermion can be represented by any 
linear combination: 



^qz 



f 1 - < 



i(f> 



(A2) 



and any choice of 4 such combinations which mutually anti- 
commute could be associated with {b x , b y , b z , c}. However, 
this difference is not physical, as all such mappings are equiv- 
alent under SU(2) gauge transformations. 

The mapping ( |A1) > does not preserve the form of the unpro- 
jected spin operators, however. Specifically, the fermioniza- 



tion Q gives 



sf = (fyn+fifit) 
si = -<(/i/u+/i/*t) 

Sf = (fi-ffit ~ flifilj 



(A3) 



while Kitaev's Majorana fermionization stipulates: 



sf 

This gives: 



ib?Ci 



ibfc; 



it 



HI 



fl + fa) (fl 



fi\, 

fa 



(A4) 



Sf = -si-i^yi + hfa 

By ox / ft ft f f 



Sf 



-Sf + (ntf + mi-l) 



(A5) 



which, after a gauge transformation to rotate the spins and 
eliminate the extra phases, differs from the spin operators 
( |A4| i by terms which vanish under projection onto the physi- 
cal Hilbert space. It is these extra terms which lead to the fact 
that the mean-field Hamiltonian d6|l does not conserve bfb x on 
x-links (and similarly for y and z) so that it is not obvious that 
the mean-field theory captures the essentials of the spin-spin 
correlations, as it is in the Majorana description. 

However, one way to view the equivalence of the two de- 
scriptions is via the wave functions that they produce after 
projection. The Majorana projector is: 



Di 



b?b!b* Ci 



(A6) 



/£ + fir) (/Jt - fir) (fl + fa) (A fa) 



Expanding the constraint in terms of Dirac fermion operators, 
we obtain 



A - -(2n it - l)(2n 4 - 1) 
= -2(n lT + n u -l) 2 + l 

Hence imposing the diagonal SU(2) constraint 



Hi 



1 = 



(A7) 



(A8) 



automatically imposes the Majorana constraint Z)j = 1. 

Therefore, if we begin with a mean-field wave-function ex- 
pressed in terms of the spinful fermions, and project onto the 
physical Hilbert space of singly occupied states, this is equiva- 
lent to studying the same mean-field wave function expressed 
in terms of Majorana fermions, and applying the projector 



( A6 1 at each site. This gives an alternative perspective on why 
the mean-field theory is exact. 
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Appendix B: Mean-field theory of the quadratic spin modef 1. Hubbard-Stratonovich decoupfing of the Kitaev model 



Here we will review the detailed derivation of the mean- 
field Hamiltonian (mi. We will first show how to derive the 
full effective action, and then present the self-consistent mean- 
field solution. 



In the Dirac fermion basis, the 3 different types of terms in 
the Hamiltonian ([TJ are: 



QX QX 



CV CV 



QZ QZ 



fftfftfiifil + filfjlfiffjt + fftfftfftfil + fi\,fjifftfif 



ft f\ 



4 

1 " 

4 . 
l r 



fftfftfiifji ~ filfjlfitfjf + fftlftPftfii + filfjifjffif 



ft ft 



fiffftfttftt + fiifjifiifii + //f/jt/jf/jt + fii.fjlfji.fii 



ft f \ 



(Bl) 



where we have used = 1 — rijj, in the last expression. 

To decouple the 4-fermi interactions using Hubbard-Stratonovich fields, we take the Lagrangian: 



8(|*i 


2 + 


l* 2 | 2 ) 




Jx 




8(16! 


2 + 


|6 2 | 2 ) 




J X 




8(|$i 


2 + 


l^l 2 ) 




•hi 




8(16! 


2 + 


|6 2 | 2 ) 




Jy 




4(|*i 


2 + 


l$2| 2 ) 


Jz 



ft,.c. 



(/;-./> + + **a (/;•.» - 4-^) + 
+ ei (44 + 44) + *e 2 (44 - 44) + h. c 



+ 'i>i (/;-/..• + fiju ) + ^2 (4/ Jt - njji ) + h.c. 
+ z6i (44 + 44) + e 2 (44 - 44) - h.c. 



(B2) 



+ $14/^ + ^4/^+^. 



4(|ei| 2 + |e 2 | 2 ) 

Jz 



©i/ it 4 + 244 



/i.e. 



where the fields 0, are to be understood as being evaluated on the link in question, and the h.c. is the hermitian conjugate 
with all spin directions reversed. We can check that this decoupling gives back the original action by integrating out the bosonic 
fields. For example, completing the square for the first line of C x gives: 



r 



8 


$1 








J X 




"T 


8 






Jx 







(f h fit + fftfa 

(fftfi-t ~ fjlfil 



J, x 



(/;• ./v + ./•;./,) I + y (/;•. /v + /;./,.) (/;•/-• + flu 
(4/it - - y (/;•./> - /;./,.) (/:•/-• - 4w) 3 ) 



2 8 



Integrating out the factors involving $1 and $ 2 gives a con- 
stant; the sum of the remaining pieces gives: 



2. Mean-field solution 



Jx 

T 



\fiffjtfjlfH- + filfjlfftfit 



(B4) 



as expected. 

Now we proceed in the usual way for mean-field theo- 
ries: namely, the fields 9 and $ have gapped amplitude fluc- 
tuations, as well as phase fluctuations. We will thus begin 
with a mean-field solution Q a ,ij{t) = A<r,ij> $<r,ij(t) = t<?,ij 
which reproduces the quadratic fermionic spectrum of the ex- 
act solution. We then consider the fate of the fluctuations of 
both gapped amplitude modes and gapless phase modes about 
mean-field. 



At mean-field level, the relevant information contained in 



Eq. ( B2 1 is that on each link there are potentially 4 bosonic 
fields: tf associated with hopping of up spins, tj, with hop- 
ping of down spins (which formally transforms in the oppo- 
site way under time reversal), and separate superconducting 
order parameters A+, for the spin up and spin down sec- 
tors. Formally, in terms of the fields of the previous section, 
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we take: 



their fluctuations according to: 



h 

A t 
A-f 
A; 



■ ($1 + i&2) on x and y links t-f 

■ ($! — i$ 2 ) on a; and y links ty 
= (i&i + 82) on y links A^ 
= (0i + i02) on x links A^ 
= (0i — i02) on x links A^ 



($1) on z links 
($2) on z links 
(i©i - 2 ) on y 

(01) on z links 

(0 2 ) on z links 



links 



2(J 



From the Lagrangian ( B2 1, the saddle-point equations are: 



+( x ,y) _ 


£ (fjlfil) 


4 x >y) 
l i 


- 4 (/jt/*t 


t w = 




A*) 
T l - 


-^-(fjlfil) 


A^ = 

t 


(fjlfil) 


A f = 


z if ft fit) 


A<*> = 


^ (fjlfil) 


Af 


= 4 (/jf/if) 


Af) = 


-j(fftM) 


Af = 


^ (fjlfil) ■ 



(B6) 



To satisfy the mean-field conditions (B6 1, we take: 



tij,l 

tij,l 

K?4 



2 *7 a; 

Tfj 

i j, 

Tf3 



on x-links 



f on y-links 



~^ij,l 

Aii,t = ^ on z -li n ks 



on x-links 



A 



ij,t 



on y-links 



.J, 
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Ay,j, = on z-links 



(B7) 



©1m 

(B5) 



which gives the mean-field Hamiltonian 



3. Theory of fluctuations about mean field 

We now turn to the fluctuations about the mean-field solu- 
tions. Since symmetry dictates that these cannot change the 
fermionic band structure, our focus will be to describe the 
bosonic degrees of freedom in this theory, and demonstrate 
that the gauge field is in a Higgsed phase with a residual Z 2 
symmetry group. 

The Hubbard-Stratonovich decoupling introduces 4 
bosonic fields: $1,2, whose saddle-point expectation values 
are associated with fermion hopping terms; and ©1,2, associ- 
ated with the spin-triplet superconductivity. We parametrize 



&2ij 

where the functions S. 



=R ( w S 'J-> ( ''' + Pij 



(B8) 



have support on x, y, and z links 
respectively, and the top (bottom) sign is taken for edges ori- 
ented from sublattice 1(2) to sublattice 2(1). 

The physical interpretation of these fields is as follows. 
$1 is associated with the spin-rotation invariant hopping 
terms familiar from spinon decompositions of the Heisenberg 
model 34 . The phase variables ay are the spatial components 
of the gauge fields associated with the constraints (J3J; fluctu- 
ations in the amplitude of this hopping term are parametrized 
by the scalar (j>. 

The remaining terms parametrize fluctuations of a con- 
densed superfluid which breaks the SU(2) gauge group down 
to Z 2 . We combine the fields associated with 0! and 2 , 
each of which is non-vanishing at mean-field either on x or 
y- links respectively, into a single pair of scalar fields p, 9 de- 
fined on all links in the lattice. Since at mean-field, 0's expec- 
tation value generates a spinful superconducting pairing, 9 is 
the phase of a charged superfluid, and hence in the condensed 
phase becomes the longitudinal component of the correspond- 
ing gauge field, p parametrizes the (gapped) fluctuations in 
this superfluid density. 

That $2, the hopping anti-symmetric in spin, is associ- 
ated with a charged superfluid is less obvious. We will show 
shortly, however, that ($ 2 ) breaks the off-diagonal generators 
of SU(2). As these are not the same as the generator broken 
by the superconducting terms, we use a new field 9 to denote 
the phase fluctuations. 

To find the residual symmetry group, we must evaluate the 
SU(2) flux through each lattice plaquette at mean-fielcR It is 
enlightening to express the fermionic degrees of freedom in 
terms of the usual BCS spinors: 



Xq 



ft,q 
fl- 



(B9) 



which transform under gauge transformations by e ta ' <T as 

X q ^e m - s Xq ■ (BIO) 

In this basis, the spin-symmetric and spin-antisymmetric hop- 
ping terms can be expressed 

**t+4.(*i) (fit f ft + fijil - f ft fit - fjifil) = 

ih+i^j) (x\xj - x]xi) 

itt-l(ij) [f! t J ft - flifil - f]t fit + fftfil) = 

itf-l(ij) {XiVzXj - x)°zXi) (BH) 
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As promised, the first term is gauge invariant under all gen- 
erators. The effect of a gauge transformation on the second 
term is to conjugate the matrix a z by e la ' a . Hence this term 
is invariant under the U(l) subgroup comprised of rotations 
about the z axis, but not under rotations by the two generators 
o x and a y . Fluctuations in 9 are therefore associated with the 

longitudinal modes of the broken generators a * . 

The remaining U(l) symmetry is broken by the supercon- 
ducting terms. As the pairing occurs here in the spin triplet 
channel, these cannot naturally be expressed in the BCS ba- 
sis; however, they are clearly charged under the residual U(l) 
symmetry fi a e lai fi a . Hence the U{1) symmetry is bro- 
ken to the Z2 subgroup fi a => ±./w, which is the residual 
gauge symmetry of the Hamiltonian. (Indeed, the spin-triplet 
superconducting terms are certainly not gauge equivalent to 
the terms associated with fct--i, guaranteeing that the SU (2) 
gauge symmetry is fully broken to Z2, rather than to a resid- 
ual U(l) as might otherwise be the case). As usual the phase 
fluctuations 8 can be absorbed, by means of a gauge transfor- 
mation, into the longitudinal modes of the broken U(\) gen- 
erator. 

As an aside: (Eq. (B8 1 reveals that the longitudinal modes 



of the broken generators are confined to x — y chains and z 
links in the lattice respectively. Since the corresponding gauge 
fluctuations are no longer purely transverse in the condensed 
phase, this means that only the residual Z2 gauge field and the 
amplitude fluctuations are free to propagate in both dimen- 
sions of the lattice. This explains, to a large degree, why the 
effect of including these bosons in the theory is so innocuous.) 

In summary, the fluctuations about mean-field are described 
by the real scalars p, p, and <fr, describing fluctuations in the 
amplitudes of the various condensed bosonic fields, and the 
SU(2) gauge field which is higgsed in a bi-adjoint represen- 
tation to a residual symmetry group Z2, which we may con- 
sider to have absorbed the remaining phase fluctuations as two 
Goldstone bosons. 

Appendix C: Mean-field theory of the gapped B phase 

Here we describe the mean-field theory in the presence of 
the 3-spin interaction which leads to the gapped topological 
B phase. We will show that the band structure discussed in 
Sect. [V]is, up to irrelevant operators, a saddle-point of an ap- 
propriate action, and thus constitutes at least a self-consistent 
mean-field solution to the fermion problem, if not a global 
minimum of the action. 

We begin by re-writing the 3-spin interaction as a sum of 
products of 6-fermion interaction terms: 



QX qV qz 
Dj D fe 



g [fiffftfilfU- f*tfftfj±fi± (CI) 



where we have used = 1 — to express S z in terms 
of down spins only. Of the possible fermion bilinears, only 

(fitfrt)>(fli.fjl)> and U)jki) (together with their ana- 
logues in the particle-particle and hole-hole channels) have 



non-vanishing expectation values at mean-field. ( (fhfkl ~ 

/fe-j-/fct) = 0). This gives us two possible ways to replace two 
of the 3 fermion bilinears by their mean-field values. First, we 
may take: 



(fyh)(fjM-(f*rfit)(fhfl) 



(C2) 



1 



which vanishes in the mean-field solution relevant to the Ki- 
taev model as the fermion bilinears are purely imaginary in 
position space. The only remaining possibility is: 



(/if/jf) (//sJ./i4-)/*4-/' s 4- (fit f ft) (fjlfkl)fklfil 

(fi^fft)(fkifft)fii.fH + (fitfft)(fji,fn)fki,fn, 

(fi-tfft)(fklfjl)filfH ~ (fi-tfft)(fftfH)fklfii- 



+(fftfft)(fkifji)fiifn - (fftfn)(fftfn)fkifii 



(C3) 



Taking < ij > to be an x-link and < jk > to be a z-link, and 



substituting in the mean-field values given in Eq. ( B7 1, this 
becomes: 



1 



fi\fk\ fklfil + fk\.fii fiifkl 

%7 \fli ~ /*•!•) yli ~ fH 



(C4) 



In light of the correspondence (lOi between our Dirac 
fermions and the Majorana basis originally used to diagonal- 
ize the problem, this is exactly the term originally proposed 
by Ref. 1 to break T and open a gap in the B phase. 

Before analyzing the resulting band structure, let us under- 
stand why we may simply replace the fermion bilinears by 
their mean-field values, as we have blithely done above. In 



fact, we can modify the Lagrangian ( B2 1 to produce just such 



a term at mean-field level. To see why this is so, we consider 
the action: 



= xlxi + X2X2 + iJ'xlxlxl + h - c - 



(C5) 



We will show that Cjr is well-approximated by the Hubbard- 
Stratonovich-like action: 

C = -|$i| 2 - |$ 2 | 2 +X2$2 

-W ($!$ 2 - xl$i - X f i* 2 ) xt + h.c. (C6) 



where Xi,2,3 are fermion bilinears. (The Lagrangian ( |B2| is of 
the general form of the quadratic terms in Eq. ( |C6[ ), albeit with 
more different scalar fields. This multiplicity of indices will 
not affect our qualitative result). The saddle-point equations 



$ 2 



*^X3($a-Xa) 
xl -iJ'xzi&i -Xi) 



(C7) 



X 



For J' = 0, the saddle-point equations specify that (f^ 
This is also the unique solution of the saddle-point equations 
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for J' 7^ (though in this case one might worry about instabil- 
ities which tend to drive $12 towards 00 if (xz) 7^ 0). Hence 
the extra term does not modify the structure of the mean-field 
equations, except inasmuch as (xi.,2) might be modified by 
the new interaction. 

As in the standard Hubbard-Stratonovich decoupling, we 
would like to integrate out $ x 2 to obtain Cj=. As the La- 
grangian (C61 is no longer quadratic in the variables &i,Xi> 



we will not be able to perform the integral exactly; rather, we 
will obtain Cjr as the lowest-order term in an expansion in J'. 
To see this, it is helpful to re-express C as: 



This means that, to lowest order in J', the effect of the 3- 
spin interaction is, exactly as originally postulated by Ref. |TJ 
to modify the band structure by adding next-nearest neighbor 
quadratic couplings. (We now also have to contend with the 
4 fermion interactions; however, when the quadratic problem 
has no Fermi surface, we do not expect these to be associated 
with instabilities of the free fermion problem and hence we 
can safely drop them without altering the qualitative nature of 
the physics.) 



£ — —IcjjjJ 2 _ |$ 2 | 2 — iJ'$!$ 2 X^ + h c + Cjr (C8) Form of the mean-field Hamiltonian with 3 spin interactions 



where = <& L — x\- I n the standard Hubbard-Stratonovich 
transformation there would be at this point no cross-terms 
coupling fermions to the scalar fields. We could therefore in- 



tegrate out the latter exactly and this prove that ( C6 1 is exactly 
equivalent to Cjr. Here we are unable to eliminate the cross- 
term &1&2X3 by further shifting the scalar fields, so that in- 
tegrating out the <f> fields will not reproduce Cjr exactly. If 
we take J' small, however, we may consider the effect of the 
cross-term perturbatively, and ask what the undesired addi- 
tions to the fermionic action will be. The exact correction is 
given by evaluating the series: 



SC 



T 



log 



i/|*i| 3 +|* 3 | 2 



n=0 



(C9) 



Terms with n odd integrate to since the action contains 
only even powers of <f>j. Hence the leading correction is of 
order J' 2 ; to linear order in J', then, we have recovered ex- 
actly the fermionic action we wanted. Since the scalar-scalar- 
fermion bilinear interaction is decidedly irrelevant (all scalars 
here are massive), we may conclude that the difference be- 
tween the action \C6\ and the true fermionic action Cjr is 
unimportant at least for the low-energy physics. 

The general form of this correction is simple to understand. 



The leading-order correction in the series ( C9 1 is proportional 
— X3X3- If we take X3 to have the form /^/ fe +, then we 



to 



have X3X3 — (xlxsY — ^ij-^fcl f° r a U r > an d a U terms in the 
series induce the same type of 'extraneous' interaction, which 
is to induce a second-neighbor 'Coulomb repulsion' term. 

We conclude that at least the low-energy structure of the 
phase we are interested in can be obtained by studying the La- 



grangian ( C6 1. We may now proceed as in Sect. B 2 obtaining 
a mean-field solution which satisfies: 



(*<) = (x\) 



(CIO) 



As noted above, the mean-field consistency conditions are 
identical to those at J' = 0; the only new feature of this saddle 
point is that it now includes quadratic terms coupling fermions 
on the same sublattice, such as: 



J($l)($2)/i/ ( 



•t 

*4 



QX QV QZ 

°i D j D k 


if r ik 


= h 


QV QX QZ 

D i D k 


if nk 


= h 


QX QZ QV 

°i °j D k 


if r ik 


= X 



Here we will derive the expression ( 50 » for the terms in- 
duced by the set of all 3-spin interactions at mean-field. There 
are three distinct 3-spin interactions that we must consider: 



(C12) 



The contributions to mean-field involve decoupling the result- 
ing 6-fermion interactions into combinations of a pair of 2- 
point functions multiplying a fermion bilinear. 

First, we show that only contributions multiplying bilin- 
ears of the form /i ff /jfe CT , fj a fka, etc., are non-vanishing. The 
mean-field eigenfunctions imply that (Sf ) = on each site. 
To show that (SfSf) = if a ^ a', we first note that if 
a — x,y and a' — z, any grouping of the resulting 4-fermion 
interaction into pairs involves one term in each pair which 
contains both a spin up and spin down fermion. Since the 
2-point functions of all terms involving spin flips are strictly 
0, these terms consequently all vanish. If a = x, a' = y, then 
we have: 



- Kfiffil + filfit)(fjffjl - fjjjt 

* ((fy f n )(hfu) + (fifn){fy*t) - h - c 



(Cll) 



(C13) 



which vanishes since the 2-point function on every link is 
purely imaginary, so that the products shown are purely real. 

The only remaining possibility is terms in which the 2- 
point functions whose mean-field expectation we take involve 
fermion operators from all 3 sites. Since all 2-point functions 
between sites i and k vanish at mean-field (this is guaranteed 
by the discrete symmetries C and T), the only posibility is 
terms which multiply fermion bilinears which couple the sites 
i and k. 

Our next task is to understand the precise form of these 
terms. For = li^, it is convenient to write Sf = fn.fl> — 

fiiH' for r n = we write Sf = // f /it - /it The re- 
sulting expressions contain couplings only between the spin- 
down fermions on sites i and k. Thus the 3-spin interaction 
does not modify the band structure of the spin-up fermions, 
which remain localized, at least at the mean-field level. 
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The quadratic couplings between the down spins induced by the 3-spin interactions can be expressed: 



QX qV qz 
D, ; dj D k 



qU qx qz 
°i °j °k 



QX QZ QU 

D i Oj O k 



^ijk;i + ^tjk;l)fkl-fi-l + C^ijk;i + ^ijk-i) fkifii + C^jfc ; ^ + ^ijfe;4.)/ij./fc-l- + (^ijk;l + ^ijfe;4.)/t|./I|. 

Ti( 7 ) ^ ft ft 



y^ijk-i 1 i]k;i)J kiJ'il ^ y^ijkii 1 ijk-i)JkiJii (-L ijk;i 1 ijk-\l J ii-l K 
(-^ijfc;t _ ^ijl^fklfu + (^ijfc;T ~~ ^ijk;^) f^ifa ~ (^ijfe;t _ ^ijk-f^ filf k l ~ (^ijfc;t 



-.(3) 



,(2) 



(4) 



(8) 



ft 



(5) 



(7) 



)fi\.fk\. 



(C14) 



with 




T (i) _ 

ijk;a 


(ftffjf) (fjcrfka) 


T (2) _ 
ijk\a 


~ (ftffjf) (fka-fj") 


T (3) 
ijk\a 


-{fltfjtKfjM 


T (4) 
ijk\a 




T (5) _ 
ijk;cr 


(fj-ffif) {fjcrfka) 


T (6) _ 


(fftMifkafjo-) 


T (7) _ 


(fjtfif)(fjafka) 


T (8) _ 


(fitfit)(fUl) 



is: 



1G 



JijJjk 
16 

JijJjk 
16 

«/ </ j A; 



16 



16 



(Af)*(^>y 

(4«>)*t«> 

(AW))* 



2iJ' 



+ (v)*A z ) ft f , a(v)*A z )* f f 



a(^)+( z ) ft ft 



with 



.J, 



16 



16 



(C18) 



(C19) 



7 7 •! f 

U IJ u jk 

16 

u •'J 

16 



t (ij)^(kj) 



For r i:) = x, we have J'SfSjS k , with ij an x-link and jfc a 
y-link. This gives the interaction: 



to) 



jk 



.Ato)/(fci) 

17 ^ 7/c 



iJ' 
+t\ x) *A 



( A (*)* A W* _ ^)* t (^ A y 4 + ( A (*)* t W* 

MM + ( — 4*^4^ 



a[ x) a\ v} 



16 



A 



to) 



ij'Jjk 



+ (f)Af* + Af)^)/i/t, 



/i|.M 

(C20) 



with 



(C15) 



where we have used t ( i k > = t (k i\ A^> k > = A^\ (Here 
we have defined A^" 6 ) = A^ ,z ' on x and z links, and — A^ 



on y links, in accordance with Eq. (B7 1 ). 

We next substitute in the mean-field values given in Eq. 
( |B7| i for t. Aon each link. We take t to be the hopping from 
sublattice 1 to sublattice 2 ( t„ — (/t^ /# 2o .)), and sim- 
ilarly for A. Here we write the induced quadratic couplings 
between two sites on sublattice 1; the couplings between sites 
on sublattice 2 are the same, but with — > — r^. 

For Tij = l\, the interaction is of the form J'SfSjS k , with 

tik) 

ij an cc-link and jk a z-link. We thus have Aj = 0, giving 
an interaction of: 



(x) _ 

t 


A J x. 

' 16 


Ai a) = 

T 


16 


■ 

t ~~ 


A J x. 

1 16 




. J u 

—I — 
16 



(C21) 



In all 3 cases, we obtain the mean-field interaction: 



±2iJ' 



fk^fii fkifii + fj^fkl fi\,fi 



t ft 
kl 



±2iJ' (fl-fki) (A- fa) 



(C22) 



2iJ' 



(x)* ,(z) ft f \{x)*,{z)* f f 

+4 s) *i z) vu,+Aftf/t/t 



(C16) 



We see that this induces a coupling only between Majorana 
modes in the dispersing band, leaving the band structure of 
the Majoranas localized on the z-links unaltered. 

Hence, the net effect of adding the 3-spin interaction, at 
mean-field level, is exactly to add the next-nearest neighbor 
couplings to the dynamical Majorana modes, while leaving 
the localized modes unchanged. 



with 



A 0*0 _ 
t _ 



16 



t {X) = 



16 



Az) Jz 
t) = —I — 

4- s 



■ (C17) 



Similarly, for = l 2 , we have J' SfSjS k , with ij a y-link 
and jk a z-link. Hence again A 



0, and the interaction 



Appendix D: Inducing Chern-Simons terms by integrating out 
fermions in the gapped B phase 

Here we will consider the 1-loop perturbative correction 
to the effective U(l) gauge field propagator due to the low- 
energy fermions in the gapped phase. We demonstrate that 
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though the Dirac point is intrinsically a property of the band 
structure of the superconductor - such that the electron bub- 
ble has both particle- particle and particle -hole contributions 
- the matrix structure about the Dirac point is such that in- 
tegrating out the low-energy fermions produces exactly the 
same Chern-Simons correction to the effective action as do- 



ing so for a normal Dirac cone. 

Since the Dirac cone is in only one of the 4 fermion bands, 
and we are interested only in the long-wavelength theory, we 
will isolate the effect of the propagator of the dispersing Ma- 
jorana band. The general form of the spin-down propagator in 
the gapped B phase is 



/ — 2m„ 



G 



Hi 



Am q 



\\-t q \ 2 



-i(A q -t q ) 2m q + iu i(A q - t q ) \ 



i(A* — t*) 2m q — iu — i(A* — t*) iu — 2m q 
2m q + iu i(A q — t q ) —2m q — iu —i{A q — t q ) 



\-^* q -t* q ) iu-2m q i(A*-t*) 



2m„ 



J 



|A 9 + ; 9 | 2 



i(A q +t q ) -iu i{A q +t q )\ 



( —iu 

-i(A*+tJ) -iu -i(A* q +t* q ) 

—iu i(A q + t q ) —iu 

\-i(A*+t*) -iu -i(A* + t*) 



i(A 



q ~ ^q) 

—iu 



(Dl) 



where we use the basis ip = 



Cql Cq2 C_ gl C_ q2 



Here we choose t q = -2J Z - J x e iq ~' h - J y e %q ~' l2 ,A q = 

J x e lq ' 11 + J y e lq l2 . In this case the bottom line is the propa- 
gator of the flat band (energies given by ±\t + A| = ±2J Z ; 
the top line is the propagator of the dispersing band, which 
captures all of the low-energy physics near the Dirac cones. It 
is easy to check that cross-terms between the two spin down 
bands vanish at 1-loop order in the fermion correction, so that 
we will drop contributions of the flat gapped band entirely. 

In the vicinity of the Dirac cone q = ( ^ , 0), at the isotropic 
point J x = J y = J z , we have 



A q - t q sa y/ZJ 



m q ~ 



(D2) 



Near this point in the Brillouin zone, then, the part of the prop- 
agator that we are interested in can be expressed as: 



Gc;q,< 



G (o) 



2 \^c-qu 



r<(sc) 

^ c;q,uj 



(D3) 

2ma z ) <g) 1 

2ma z ) ® a x 



with — ( 1 <Jy a x y In addition to the usual term 
(G^ q>u ), the fermion propagator contains an anomalous term 

(sc) 

(Gc-.q'u) due to the presence of superconductivity. The 2x2 
matrix structure of both of these terms is, however, the same. 

In this long-wavelength limit, the interaction between 
fermions and the gauge field is 



(D4) 



where 7 M = a „ ® 1, and the last term occurs due to normal 
ordering. (Here it should be understood that the sum encom- 
passes only half the Brillouin zone). The 1-loop correction to 
the gauge field effective action induced by the fermion terms 
is therefore: 



J \2ir) 



lvG c - 



q-\-p,uj-\-Q,\ 



(D5) 



Using the expression (|D3[), we find that traces of the cross- 



terms between G^\^ and G^q'^ vanish, leaving: 



(sc) 



(D6) 



d 3 p 
(2nf 



Tr 



where c]w is the effective action induced by the usual 2 + 1 
dimensional Dirac cone (appearing here with a multiplicative 
factor of 2 since we have counted both terms of the form 
f qi f q i and f_ q if-qi, effectively counting the contribution of 
both Dirac cones). The second contribution, due to the su- 
perconducting terms, also has precisely the same form as the 
first, since G^ has the same 2x2 structure as G^°\ The 
factor of \ (due to the \ in G^ relative to its usual value) is 
exactly cancelled by the factor of 4 from these contributions. 
This gives exactly the 1-loop correction expected from a sin- 
gle Dirac cone in QED, albeit with a mass of 2m rather than 
m. 
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